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In a superconductor, the ratio of the carrier density, n, to their effective mass, 
m* , is a fundamental property directly reflecting the length scale of the superfluid 
flow, the London penetration depth, Xl- In two dimensional systems, this ratio n/m* 
(~ 1/Af,) determines the effective Fermi temperature, 2>. We report a sharp peak in 
the ^-dependence of Xl at zero temperature in clean samples of BaFe2(Asi_ I P I )2 at 
the optimum composition x = 0.30, where the superconducting transition temperature 
T c reaches a maximum of 30 K. This structure may arise from quantum fluctuations 
associated with a quantum critical point (QCP). The ratio of T c /Tf at x = 0.30 is 
enhanced, implying a possible crossover towards the Bose-Einstein condensate limit 
driven by quantum criticality. 
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In two families of high temperature superconductors, 
cuprates and iron-pnictides, superconductivity emerges 
in close proximity to an antiferromagnetically ordered 
state, and the critical temperature T c has a dome shaped 
dependence on doping or pressure [H-Q- What happens 
inside this superconducting dome is still a matter of de- 
bate In particular, elucidating whether a quantum 
critical point (QCP) is hidden inside it (Figs.[TJ\ and B) 
may be key to understanding high-T c superconductivity 
[3, Hj]. A QCP marks the position of a quantum phase 
transition (QPT), a zero temperature phase transition 
driven by quantum fluctuations Q. 

The London penetration depth Ar, is a property that 
may be measured at low temperature in the supercon- 
ducting state to probe the electronic structure of the 
material, and look for signatures of a QCP. The absolute 
value of Xl in the zero-temperature limit immediately 
gives the superfluid density A^ 2 (0) = /ioe 2 JY rii/m* , 
which is a direct probe of the superconducting state; here 
m* and rij are the effective mass and concentration of the 
superconducting carriers in band i, respectively 0. Mea- 
surements on high-quality crystals are necessary because 
impurities and inhomogeneity may otherwise wipe out 
the signatures of the QPT. Another advantage of this 
approach is that it does not require the application of a 
strong magnetic field, which may induce a different QCP 
or shift the zero-field QCP [|. 

BaFe2(Asi_ 2; P 2 ;)2 is a particularly suitable system for 
penetration depth measurements as, in contrast to most 



other Fe-based superconductors, very clean [l(| and ho- 
mogeneous crystals of the whole composition series can 
be grown 11|. In this system, the isovalent substitu- 



tion of P for As in the parent compound BaFe2As2 of- 
fers an elegant wa y t o suppress magnetism and induce 
superconductivity [111 ]. Non- Fermi liquid properties are 
apparent in the normal state above the superconducting 
dome (Fig. [2^1 [11 E2] and de Haas- van Alphen (dHvA) 
oscillations [10] have been observed over a wide x range 
including the superconducting compositions, giving de- 
tailed information on the electronic structure. Because 
P and As are isoelectric, the system remains compen- 
sated for all values of x (i.e., volumes of the electron and 
hole Fermi surfaces are equal). 

As discussed in Ref.llOl the normal-state electronic 
structure of BaFe 2 (Asi_ a ;P ;E )2 determined by dHvA ex- 
periments is significantly modified from that predicted 
by conventional density functional theory (DFT) band 
structure calculations. Figure [3JA shows the composi- 
tion evolution of the effective mass m* normalized by 
the free electron mass m and the Fermi temperature 
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-Ah for x > 0.4, determined from 



the dHvA oscillations corresponding to the extremal or- 
bits on the outer electron Fermi surface (ft and ft orbits 
in Fig.[2j3). Here e_F is the Fermi energy and Ak is the 
cross sectional area of the orbit. In contrast to the negli- 
gible x-dependence expected from the DFT calculations, 
a critical-like increase in m* accompanied by a strong re- 
duction of Tp is observed as the system is tuned towards 
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FIG. 1. Generic temperature vs. nonthermal control param- 
eter phase diagram of iron-based superconductors illustrating 
two cases. (A) Quantum criticality is avoided by the transi- 
tion to the superconducting state. There is only one super- 
conducting phase. (B) A QCP lies beneath the superconduct- 
ing dome. The QCP separates two distinct superconducting 
phases (SCI and SC2). In the case of A, non- Fermi liquid be- 
havior may appear above the dome if there is a QCP located 
along the axis of another control parameter that is indepen- 
dent of the control parameter shown on the abscissa. In the 
case of B, non- Fermi liquid behavior appears due to the QCP 
inside the dome. 



the optimal composition from the overdoped side. 

For a reliable determination of the absolute value of 
Al(0) in small single crystals, we adopted three different 
methods @. The first is the lower- T c superconducting 
film coating method [l3l - [l5l |. in which Ai(0) is deter- 
mined from the frequency shift of a high precision tunnel 
diode oscillator [T(| (resonant frequency of / ~ 13 MHz) 
containing the BaFe2(Asi_ 2; P a ;)2 crystal coated with an 
aluminum film (T c — 1.2 K) of known thickness and pen- 
etration depth. 

The second is the microwave cavity perturbation tech- 
nique, in which A^(0) is determined from the measure- 
ments of surface impedance, Z s = R s + iX s , by using 
a superconducting resonator (/ ~ 28 GHz) and a rutile 
cavity resonator (/ ~ 5 GHz), both of which have a very 
high quality factor Q ~ lfj 6 J8[. In all crystals, the resid- 
ual surface resistance R s (0) at T — >• OK, which we deter- 
mined by withdrawing the crystal from the rutile cavity 
at low temperature, is less than 0.3% of R s just above T c . 
This negligible residual R s (0) indicates almost perfect 
Meissner screening without any non-superconducting re- 
gions. In the superconducting state well below T c , \l(T) 
is obtained from the surface reactance via the relation 
X S (T) = /LtowAi(T). The absolute value of X s is de- 
termined from Z s and dc-resistivity pd c (measured sepa- 
rately by a conventional four contact technique) by the 
relation R s = X s = \J fiQUjpdc/2 which holds in the nor- 
mal state [8j. 

The third method uses the temperature dependent 
changes SXl{T) — Al(T) — Al(0), measured by the tunnel 
diode oscillator down to ~ 80 mK (Fig. For all sam- 
ples measured, covering a wide range of concentrations 
0.26 < x < 0.64, a quasi-T-linear variation of SXl{T) is 
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FIG. 2. (A) Phase diagram of BaFe 2 (Asi_ ;r P ;c )2. The tran- 
sition to the SDW ground state at Tm coincides or is pre- 
ceded by the structural transition at T s . With increasing 
x, Tn decreases and goes to zero continuously at x = 0.30. 
The superconducting dome extends over a composition range 
0.2 < x < 0.7, with maximum T c = 30 K at x = 0.30. The 
red shaded region at around x = 0.30 represents the region 
where the exponent n of the temperature dependence of the 
resistivity, pdciT) — p(0) + aT n , is close to unity, which is 
a hallmark of a non- Fermi liquid (non-FL). The composition 
dependence of the effective Fermi temperature 2> and renor- 
malized mass m* /mo determined by dHvA oscillations fiol ] 
arising from the /3 orbits (shown in B) are also plotted. (B) 
Fermi surface of BaFe2(Asi_ :E P a ;)2 with x = 0.3 and 0.7 from 
the band-structure calculation using density functional the- 
ory as implemented in the WIEN2K code [13]. The Fermi 
surface consists of five quasi-cylindrical pockets, three hole 
pockets at the center of the Brillouin zone, and two electron 
pockets centered at its corners. The shading represents the 
in-plane Fermi velocity vf. The flat parts of the outer elec- 
tron sheets have high vf values. The lines represent the ex- 
tremal /3-orbits. (C) Composition evolution of the square of 
the London penetration depth A^(0) in the zero-temperature 
limit determined by three different methods: Aluminum coat- 
ing method (black diamonds) , microwave cavity perturbation 
technique (blue circles), and the low-temperature slope of 
the change of the penetration depth with temperature (red 
squares, right hand scale) shown in Fig. 3. Different points 
for the same x correspond to different crystals from the same 
batch. Error bars shown for x — 0.64 represent typical exper- 
imental errors for all x, which involve the measurement errors 
and typical systematic uncertainties including the sample 
dependence. 



observed. This important result indicates that the pres- 
ence of line nodes in the superconducting gap [ll| is a 
robust feature of this P-substituted system. This robust- 
ness is consistent with the nodes being on the electron 
sheets 17] rather than the hole sheets, as the electron 
sheets change relatively little with x whereas the shape of 
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FIG. 3. Relative change of the penetration depth, SXl(T) = 
Al(T) — Al(0), at low temperatures plotted against T/T c for 
different compositions from x = 0.27 (A) to 0.64 (F). For 
x = 0.30, data for two samples are shown, one of which (#2) 
is shifted vertically for clarity. 



the hole sheets changes substantially (Fig. 03). A notable 
feature of the T-linear penetration depth is that the rela- 
tive slope S\ L /d(T/T c ) is steepest for x = 0.30 (Fig. [3^). 
In general, this slope is determined by the Fermi veloc- 
ity and the k dependence of the superconducting gap 
close to the node. Making the reasonable assumption 
that the gap structure evolves weakly across the phase 
diagram, the x dependence of 5Xl/<1(T/T c ) will mirror 
that of A L (0) 0. 

Figure [2C shows the composition dependence of the 
squared in-plane London penetration length A 2 (0) in the 
zero-temperature limit. Although different techniques 
can involve systematic errors [l5[ , all three methods give 
very similar x dependencies. The most notable feature is 
the sharp peak in A|(0) at x = 0.30, at about the same 
composition level where T c is maximal. The striking en- 
hancement of A^O) is observed on approaching x = 0.30 
from either side, and has been seen in multiple samples 
using different techniques. This reproducibility, com- 
bined with the above mentioned low R s (0), sharp super- 
conducting transitions, and large heat capacity anoma- 
lies at all values of x close to x = 0.30 [|[, shows that the 
enhancement is not an experimental artifact associated 
with poor screening caused by non-bulk superconductiv- 
ity. We attribute the peak in A|(0) to the existence of a 
QCP at x = 0.30. 

This result contrasts with the behavior found 
for the electron-doped iron-based superconductors, 
Ba(Fei_ x Co x )2As2, where a shallow minimum of Al(0) 
at the optimum doping and a continuous increase on the 
underdoped side have been reported EH, HH . This dif- 
ference from the present case may be related to a greater 
degree of electronic disorder in the Fe layer caused by 
the Co doping [IS], which may smear out the singular- 
ity. The difference in the superconducting gap structure 
[2| as well as the additional of charge carriers from the 
Co doping, may also be a source of differences in the 



x-dependence of Al(0). 

In cuprates, a QCP associated with the pseudogap for- 
mation has been suggested to occur at the hole concen- 
tration p ~ 0.19 inside the superconducting dome. We 
note, however, that there does not appear to be any evi- 
dence of mass divergence at this purported QCP and at 
this doping a broad minimum of A^(0) was reported [22j . 
An enhancement in A 2 (0) has been observed at p ~ 1/8 
[23j . but this is accompanied by a reduction of T c , which 
is again different to the present case where the peak in 
A|(0) coincides with the maximum T c . 

Our results may have general implications for the be- 
havior of A^O) in strongly correlated superconducting 
systems. How strong electron correlations influence the 
condensed electron pairs in superconductors has been a 
long-standing issue |24l - t2a |. In fact, it has been pointed 
out that in an ordinary one-component Galilean invari- 
ant Fermi liquid, electron correlation effects do not cause 
the renormalization of Al in the superconducting state 
j24j. However, experimentally A^O) does appear to be 
enhanced in heavy-fermion superconductors, which con- 
tain interacting conduction electrons and local moments 
[25|,[26j]. The present results in BaFe2(Asi_ a; P :I ;)2 sup- 
port this and suggest that in sufficiently clean systems 
electron correlation effects can lead to a striking renor- 
malization of A|(0). 

We now discuss the consequences of a QPT inside the 
superconducting dome. Such a QPT implies that the 
non-Fermi liquid behavior indicated by the red region in 
Fig.[5K is most likely associated with a finite tempera- 
ture quantum critical region linked to the QCP. More- 
over, this transition immediately indicates two distinct 
superconducting ground states. In our system, the ro- 
bust T-linear behavior of 6 A l (T) on both sides of the pur- 
ported QCP at x = 0.30 argues against a drastic change 
in the superconducting gap structure 0, |3{|. The fact 
that the zero-temperature extrapolation of the antiferro- 
magnetic transition Tv(x) into the dome [l2| coincides 
with the location of the QCP (Fig.[2j\) may indicate that 
the QCP separates a pure superconducting phase on the 
right and a superconducting phase coexisting with spin- 
density-wave (SDW) order on the left (Fig.Qj3). 

To place BaFe2(Asi_ 2 ,P a; )2 in the context of other su- 
perconductors, Fig.0] plots T c as a function of the ef- 
fective Fermi temperature Tp for several types of com- 
pounds (Uemura plot); the red symbols correspond to 
various values of x for BaFe2(Asi„ x P2;)2 studied here, 
and the others are obtained from /iSR measurements 
reported previously [27j . Because the relevant Fermi 
surface sheets are nearly cylindrical, Tp for 2D sys- 
tems may be estimated directly from Al(0) via the re- 



lation, Tp = 



(g_7r)w2D 



( U 2 TC 



A L 2 (0), where ri2D 



is the carrier concentration within the superconduct- 
ing planes and d is the interlayer spacing; Tp = 
(h 2 /2)(35T^) 2/3 n 2/3 /k B m* for three dimensional (3D) 
systems (27]]. The dashed line corresponds to the Bose- 
Einstein condensation (BEC) temperature for an ideal 
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quasi-2D system, this value of Tg provides an estimate of 
the maximum condensate temperature. The evolution of 
T c with Tf in the present system is in sharp contrast to 
that in cuprates, in which T c is roughly scaled by Tp ■ The 
inset of Fig. reffig4 depicts the ir-composition dependence 
of Tc normalized by Fermi (or BEC) temperature, T c /Tf 
(T C /T B ). In the large composition region (x > 0.6), 
T c /Tp is very small, comparable to that of the conven- 
tional superconductor Nb. As x is decreased, T c /Tp in- 
creases rapidly, and then decreases in the SDW region 
after reaching the maximum at the QCP (x = 0.30). 
Notably, the magnitude of T C /T B (« 0.30) at the QCP 
exceeds that of cuprates and reaches almost 40% of the 
value of superfluid 4 He. 



FIG. 4. Uemura-plot. T c is plotted as a function of effective 
Fermi temperature Tf evaluated from 1 /\ 2 L (0) for various su- 
perconductors (w2p/(m*/mo) for 2D and 1.52n/(m*/mo) for 
3D systems) [23] . We used an average of the Al-coating and 
microwave data for BaFe2 (Asi- x P x )2- The data for x > 0.30 
(red circles) and for x < 0.30 (red squares) bridge a gap 
between the conventional superconductors such as Nb and 
cuprate high-T c superconductors such as (La,Sr)2Cu04 (214), 
YBa2Cu 3 7 -5 (123), and Bi2Sr 2 Ca2Cu 3 H (2223). x = 0.30 
represents the data at the QCP. The dashed line is the BEC 
temperature for the ideal 3D boson gas. Inset: composition 
dependence of T c normalized by the Fermi temperature (left 
axis) or BEC temperature (right axis). Green and light blue 
dashed lines mark the T c /Tf values for underdoped cuprates 
123 and for the conventional superconductor Nb. Brown ar- 
row represents T c /Tb = 0.7 for superfluid 4 He. 



The fact that T c /Tp becomes largest at the QCP in- 
dicates that the strongest pairing interaction is achieved 
at the QCP, implying that high-T c superconductivity is 
driven by the QCP. In a multiband system, we need to 
introduce the effective Fermi energy tp for each band, 
which is defined for electron (hole) bands as the energy 
of the highest occupied state relative to the bottom (top) 
of the band. Because the outer electron sheet with the 
highest Fermi velocity has the largest tp and hence the 
largest contribution to A^ 2 (0), the magnitude of T c /Tp 
in the other sheets are expected to be even larger. These 
results lead us to consider that in terms of TJTp the 
system is closer to the BCS-BEC crossover [28 30] than 
the cuprates. 
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I. SAMPLE CHARACTERIZATION 

High-quality single crystals of BaFe2(Asi_ a; P 2 ,)2 were 
grown by the self-flux method [11]. For the crystal struc- 
ture analysis, X-ray diffractometry was carried out using four- 
circle diffractometer (MXC X , MacScience). We find that the 
composition dependence of the lattice parameters a, c, and 
the pnictogen height hp n obey Vegard's law, which allows us 
to determine the P-concentration x within the error of ±0.006. 
In this study, BaFe2 (Asi-aPa,^ single crystals in a wide range 
of P-concentrations (0.23 < x < 0.64) were grown for the 
measurements of the London penetration depth. Typical size 
of the crystals is ~ 200 x 200 x 10 pm 3 . The crystals were 
cleaved just before the measurements to obtain the clean and 
flat surfaces. Figure lSll shows the temperature dependence of 
the microwave resistivity pi = 2R^/po<-^ at 4.9 GHz for the 
crystals used in this study, which coincides with the dc resis- 
tivity pdc in the normal state. For x — 0.27, a kink anomaly 
associated with the structural or SDW transition is observed 
at 47K, which is consistent with the previous studies [11]. 
For x = 0.30 at the optimal P-concentration, pi shows a T- 
linear dependence, which is the hallmark of non-Fermi liquid 
behavior. For x — 0.64, pi shows a T 2 -dependence, indicating 
a recovery of the Fermi- liquid behavior in the high concentra- 
tion regime. All these pi (T) behaviors are consistent with the 
previous pdc measurements [11]. Figure [S2l shows the normal- 
ized ac susceptibility measured by the tunnel diode oscilla- 
tor. A very sharp superconducting transition is observed for 
x > 0.30, demonstrating the high quality and homogeneity. 
The almost perfect superconducting volume fraction in the 
sample of x = 0.30 is supported by the large heat capacity 
anomaly (inset of Fig. lS2|) . The superconducting transition 
becomes broader for x < 0.30, because the steep slope dT c /dx 
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FIG. S2. Normalized ac susceptibility of BaFe2(Asi- :r P; C )2 
with different P concentrations x (0.26 < x < 0.64) measured 
from the frequency shift of the tunnel diode oscillator. The 
inset shows the specific heat C with the normal state contribu- 
tion subtracted for the same x = 0.30 sample as was used for 
the temperature dependent penetration depth measurements 
(Fig. 3). The normal state contribution was estimated by ex- 
trapolating C/T measured above T c . The measurements were 
made using a modulated temperature method [32] , which has 
high resolution but poor absolute accuracy for small samples. 
The absolute value of AC was therefore estimated by normal- 
ising the total specific heat to that measured on a similarly 
doped sample by a dc technique [18] . The uncertainty in AC 
is ~ ±20% which is determined mostly by uncertainties in the 
addenda contribution. 



in the phase diagram (Fig. 2A) implies that T c is very sensitive 
to small variations of x in the low concentration regime. The 
bulk T c is defined by the temperature at which superfluid den- 
sity becomes zero extrapolated from low temperatures. We 
note that this T c coincides with the temperature at which spe- 
cific heat and thermal expansion coefficients jump [31] and dc 
resistivity vanishes. 
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FIG. SI. Low-temperature microwave resistivity pi of 
BaFe2(Asi_ I P a: )2 at 4.9 GHz for various P-concentrations. 
The arrow indicates a kink anomaly observed for x = 0.27, 
which is attributed to the structural or SDW transition. The 
lines are guides for the eye. 



II. EXPERIMENTAL TECHNIQUES 
A. Tunnel diode oscillator technique 

The tunnel diode oscillator (TDO) technique is able to 
extract the relative change of the penetration depth with 
temperature, SX(T) = A(T) — A(0). The penetration depth 
measurements down to very low temperatures have been per- 
formed by using a high resolution radio frequency susceptome- 
ter based on a self-resonant tunnel diode circuit operating at 
13 MHz [33], which is mounted on a 3 He- 4 He dilution refrig- 
erator (base temperature ~ 80 mK) or 3 He refrigerator (base 
temperature ~ 300 mK). 

The TDO operates with an extremely small ac probe field 
(Hac < lOmOe) so that the sample is always in the Meiss- 
ner state. To measure the in-plane penetration depth, H ac is 
applied along the crystal c axis, which generates the super- 
currents in the plane. The change in the resonant frequency 
5f is proportional to the change in the penetration depth 8\, 
8f = G5X. The calibration factor G is determined from the 
geometry of the sample, and the total perturbation to the 
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resonant frequency due to the sample, found by withdrawing 
the sample from the coil at low temperature [34]. The sample 
is mounted on a sapphire rod, the other end of which is glued 
to a copper block on which a RuC>2 or Cernox thermometer 
is mounted. The sample and rod are placed inside a solenoid 
which forms part of the resonant tank circuit. 



B. Microwave cavity perturbation technique 

The penetration depth can also be determined by the 
measurements of the microwave surface impedance, Z a = 
R a + iX s , where R a is the surface resistance and X 3 is the 
surface reactance. We have measured Z a by high-resolution 
microwave surface-impedance probes that use cavity per- 
turbation of superconducting and dielectric resonators. A 
cylindrical superconducting cavity resonator with operating 
frequency of 27.8 GHz in the TEon mode is made of Pb- 
plated oxygen free high conductivity (OFHC) copper [35,36]. 
A cylindrical dielectric cavity resonator with operating fre- 
quency of 4.9 GHz in the TEon mode is constructed from a 
high-permittivity single crystal of rutile (Ti02) with high di- 
electric constant (e w 120) [37]. Both resonators have quality 
factors in excess of 10 6 . 

The superconducting cavity resonator is immersed in su- 
perfluid 4 He at 1.6 K and the temperature of the whole res- 
onator is stabilized within ±0.5 mK. Figure lS3K shows the 
schematic of the dielectric cavity resonator. The magnetic 
flux is mainly confined in the rutile, which reduces the con- 
duction current losses in the walls of the resonator enclo- 
sure. Figurc [S3B illustrates the magnetic field structure of 
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FIG. S3. (A) Schematic of the rutile cavity resonator. The 
rutile resonator is placed inside the OFHC copper enclosure. 

(B) Left panel: an illustration of the magnetic-field distribu- 
tion of the TEon mode. Right panel: z dependence of |//| 2 . 

(C) -(E) The residual surface resistance at 4.2 K normalized 
by the value just above T c at z = is plotted as a function of 
z for the samples of x = 0.27 (C), x = 0.30 (D), and x = 0.56 
(E). 



the TEon mode. The temperature of the rutile is stabilized 
within ±0.5 mK. The crystal was mounted on a sapphire hot 
finger and placed at the antinode of the microwave magnetic 
field (|| c axis) so that the shielding currents 1^ are ex- 
cited in the ab plane (Fig. lS3B ). The inverse of quality factor 
1/Q and the shift in the resonance frequency are proportional 
to the surface resistance R s and the change in the surface re- 
actance SX a = X S (T) — X s (0), respectively. The sample stage 
is movable along the vertical axis z of the sapphire rod from 
the center of the cavity (2 = 0). By withdrawing the sample 
completely from the resonator, the background microwave ab- 
sorption of the resonator can be measured directly, allowing 
the absolute microwave absorption of the sample to be deter- 
mined in situ. Figures [S3C -E show the z-dependence of the 
residual surface resistance R s at 4.2 K normalized by the value 
just above T c at z = for three samples with different compo- 
sitions. The residual surface resistance at 4.2 K in all crystals 
is less than 0.3% of R s at T c , demonstrating nearly perfect 
Meissner state with negligibly small non-superconducting re- 
gion. We note that even in the low-concentration samples 
with broader transitions the residual R s value is comparable 
to that in optimal composition (Fig. lS3C ). This ensures that 
Al (0) values in these samples give correct estimate of intrinsic 
superfluid density. 



III. DETERMINATION OF THE ABSOLUTE 
VALUE OF THE ZERO TEMPERATURE 
LONDON PENETRATION DEPTH 

A. Al-coating method 

To determine the absolute value of the penetration depth, 
we used the lower- T c superconducting film coating method 
[38]. We determined A(0) from the the frequency shift of 
the TDO containing the BaFe2(Asi_ a; P ;E )2 crystal whose sur- 
faces are coated with a thin film of aluminum (Al) having a 
lower critical temperature (T^ 1 = 1.2 K) and a known value 
of the penetration depth Aai(0) = 55 nm [13,14]. FigurelSlR 
shows the temperature dependence of the penetration depth 
for x — 0.3 sample coated with the Al film of thickness of 
~ 100 nm. Above T C A1 fFig.lSIBl, the normal-state skin depth 
of Al (<Sai ~ 75 fim for p A1 — 10 /ificm at 13 MHz) is much 
larger than the thickness of the Al film. As a result, the 
effective penetration depth into both the Al film and coated 
superconductor, A e ff(T), is almost identical to the penetration 
depth A(T) before coating Al on the sample. On the other 
hand, when the Al film is superconducting (Fig. lS4C ). Al acts 
together with the coated superconductor to screen the mag- 
netic fields. Then, the effective magnetic penetration depth 
Acff(T) for T < T C A1 can be given by 



Aeff(T) = Aai(T) 



A(T) + Aai(T) tanh ■ 



l(T) 



Aai(T) ± A(T)tanh 



(SI) 



Aai(t) 



where t is the thickness of the Al film and A(T) is the pene- 
tration depth of the coated superconductor and Aai(T") is the 
penetration depth of the Al film [14]. Thus we can extract 
the absolute penetration depth of the coated superconductor. 
The overall penetration depth change below the Al transition, 
L = 5\ c ff(T) — A c ff(T) — A e ff(T m i n ), is used for the calcula- 
tion of A(0), which gives A(0) = 330 nm for x — 0.3. The 
typical error bar for Al(0) measurements by this technique is 
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FIG. S5. Surface resistance R a and reactance X a as a function 
of the temperature at 4.9 GHz for x — 0.56. In the normal 
state, R s = X a is observed. Inset shows the temperature 
dependence of the absolute penetration depth. Solid line rep- 
resents the TDO data measured down to 80 mK on the same 
single crystal. 



FIG. S4. (A) The temperature dependence of the effective 
penetration depth for the optimally doped BaFe2(Aso.7Po.3)2 
sample coated with Al film. The inset shows the low- 
temperature behavior. The overall variation of the penetra- 
tion depth below the Al transition, denoted as L, is used for 
the calculation of A(0). (B and C) Schematic of the field pen- 
etration in Al-coated sample at T > T C A1 (B) and T < T C A1 
(C). 



trapolated down to zero temperature by using the TDO data 
measured down to 80 mK, which gives A(0) = 135 nm. The 
typical error bar for Al(0) measurements by this technique is 
±15%. The obtained results are summarized in Fig. 2C (blue 
circles) . 



C. Temperature dependent changes SXl(T) 



±15%. The obtained results for the other concentrations x 
are plotted in Fig. 2C (black diamonds). 



B. Surface impedance measurements 

Figure lS5l shows the temperature dependence of the surface 
resistance R a and reactance X s for the x = 0.56 sample at 
4.9 GHz. The low residual R a value in the low-temperature 
limit again demonstrates the high quality of the crystal. The 
transition in microwave R a is broader than that in dc re- 
sistivity data since the applied GHz microwave excites ad- 
ditional quasiparticles particularly just below T c . In the 
superconducting state well below T c , A(T) is obtained via 
the relation X a (T) — /iowA(T). In the normal state, at 
the frequencies lot <C 1, the surface resistance is equal to 
the surface reactance in the so-called Hagen-Rubens limit, 
R a — X a = *J noLopdc/2. As shown in Fig. IS5I the relation 
R a = X a holds in the normal state, indicating the Hagen- 
Rubens limit in the present frequency region. The absolute 
value of X a is thus determined by the dc-resistivity pdc in 
the normal state. The inset of Fig. lS5l shows the temperature 
dependence of the absolute value of the penetration depth 
for the x = 0.56 sample, demonstrating a good agreement 
with the TDO data shifted by a constant value. A(T) was ex- 



The quasi- T- linear variation of 5\ L (T) = \ L (T) - Xl(0) 
is observed in a wide range of x, indicating line nodes in the 
gap. This temperature linear slope is related to Al(0) and 
the superconducting energy gap magnitude A. For a d-wave 
nodal superconductor it is found that [13] 



S\ L (T) _ ln2 
Ai(0) ~ A 



k B T. 



(S2) 



If we assume this formula and a constant ratio A/fcsT c 
(which in the simplest case corresponds to 2.14 for d-wave 
and 1.76 for s-wave), then A|(0) can be estimated from the 
measured dX/dT and T c . For example, when we use the s- 
wave ratio for x — 0.3, dX/d(T/T c ) = 125 ± 10 nm, giving 
Al = 320±30 nm. This compares to 330 nm obtained from the 
Al coating method and 300 nm from the surface impedance 
methods. The red squares in Fig. 2C are plotted in a scaled 
right axis, such that A L (0) estimated by these simple assump- 
tions (with A/k B T c = 1.76) from (T c ■ dX/dT) 2 can be com- 
pared with the values in the left axis obtained directly by the 
other methods. Although the relation between dX/d(T/T c ) 
and Ai,(0) in general depends on the number of nodes, the 
rate at which the gap opens near the node (dA/dk) and also 
the relation between A and T c , the correspondence with this 
simple estimate is striking. In any case, as long as these fac- 
tors do not change sharply at the QCP it the observed strong 
increase in dX/d(T/T c ) close to the QCP is strong evidence 
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supporting the other measurements which indicate a sharp 
peak in Al(0) at this composition. 



IV. MICROSCOPIC COEXISTENCE OF 
SUPERCONDUCTIVITY AND SDW ORDER IN 
THE UNDERDOPED REGIME 

There has been much debate as to whether superconduc- 
tivity coexists microscopically with other ordering not only in 
pnictides but also in cuprates, heavy fermions and other ex- 
otic superconductors, because of the difficulty to distinguish 
this from mesoscopic phase separation. The present results 
that the QCP separates two distinct superconducting phases 
strongly suggest that in the underdoped region, superconduc- 
tivity and SDW coexist on a microscopic level, but compete 
for the same electrons in the underdoped region. This com- 
petition is evidenced by the overall larger Al(O) values in the 
SDW side of the QCP than the other side (Fig. 2C), corre- 
sponding to the smaller volume of Fermi surface due to partial 
SDW gapping. The microscopic coexistence is also supported 
by the enhancement of A 2 (0) on approaching the QCP from 
the SDW side, which is not expected in the case of phase 
separation [14,39]. 



V. MICROSCOPIC INTERPRETATION OF THE 
MAGNETIC PENETRATION DEPTH 

The magnetic penetration depth A is defined via the pene- 
tration of field into the surface of a superconductor 



A 



Bo jo 



Bdr. 



In a clean local, BCS superconductor, at T = 0, the x com- 
ponent of A can be evaluated by the following integral (in SI 
units) of the Fermi velocity v over the complete Fermi surface 
S [40] 



A" 2 (0) 



fi e' 



[ < 
J M 



dS. 



(S3) 



47T 3 fi, 

Here v x is the x component of the Fermi velocity. This equa- 
tion is closely related to that for the normal state conductivity, 



where r is the scattering time which in general depends on 
k. In the free electron, isotropic scattering approximation, 
these two equations can be reduced to A _2 (0) = none 2 jm and 
a = ne 2 r/m, where m is the mass of the electron and n is 
the density of electrons, and so it is common in the literature 
to abbreviate these integrals to the quantity n/m* , where 



n 

m* 



fl ds - 

\v\ 



-dS 



(S4) 



For a multiband material, a separate n/m* can be evaluated 
for each sheet of Fermi surface and the contributions added 
to give the total A^ 2 (0). The 'mass' m* here is closely related 
to the cyclotron mass m* which is measured in de Haas-van 
Alphen effect experiments 



h 2 dA 

where A is the Fermi surface (extremal) cross-sectional area 
and £ is the quasiparticle energy. Both masses are renor- 
malized by the same factor when the Fermi velocity is renor- 
malized by many-body effects. This assumes that 'back-flow' 
cancellation of the Fermi liquid renormalization of A -2 [24,25] 
does not occur in real superconductors which are not trans- 
lationally (Callahan) invariant. This is supported by experi- 
ment [26]. n is proportional to the volume of the Fermi surface 
and according to Luttinger's theorem [41] it is independent of 
the strength of the correlations. However, in iron-based su- 
perconductors where there are both electron and hole parts to 
the Fermi surface it should be noted that although the total 
volume is independent of the interactions, n in the above for- 
mula may change. This is because although the total volume 
takes into account the sign of the carriers (it is zero for a com- 
pensated metal such as BaFe2(Asi_ I P I )2 where the electron 
and hole pockets have equal volume), the sign of v does not 
enter Eq. ([S3} and so the electron and hole pocket contribu- 
tions add. The dHvA results [10] suggest that although the 
volumes of the electron and hole pockets shrink as x decreases 
from 1 towards the optimum composition, the main effect on 
n/m* near the QCP comes from the increase in m* . It should 
also be noted that A -2 may be reduced from the value derived 
from Eq. (|S3[) by phase fluctuations [42], non-locality or dis- 
order. 
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